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In a previous report, we introduced a new fermionic variational wavefunction, suitable for interact- 
ing multi-species systems and sustaining superfluidity. This wavefunction contains a new quantum 
index of "ghostly" character, which does not appear in any observable quantity. Here we intro- 
duce a spin triplet version of this wavefunction, with parallel spin pairs only. Moreover, we present 
numerical solutions for the spin singlet superconducting state previously introduced. 

o 

In a previous report []| , to be referred to as (I) , we introduced a new variational wavefunction, suitable for interacting 

^ . multi-species systems and sustaining superfluidity. This a generalization of the Bardeen-Cooper-Schrieffer (BCS) 

'""' ^ wavefunction |2| I^bcs) — Wki'^k + Vk cj, >,c_j, i )|0). The creation/annihilation operators cj, ^/ck,(T describe fermions 

lO ■ with momentum k and spin a, and |0) is the vacuum state. 

' We introduce a spin triplet version of this wavefunction, which corresponds to the "equal spin pairing" (ESP) case 

r^ . with parallel pair spins only. In principle, the ESP state does not yield the lowest ground state ^3ij for the single 

O ' species case, which may apply here as well. 

^ . Let the usual fermionic operators be cl/cx with x = {z, k, a}, where i denotes the fermion species/flavor. 

$-H ' We introduce a new quantum index, which is related to the internal symmetry of the quantum state. It serves 

^-1, to enumerate the otherwise "entangled" components of the quantum state as a function of both momentum and 

^ ■ spin. In this way, the treatment of the coherence factors (coefficients entering the wavefunction and then everywhere 
else in the theory) is greatly facilitated. Thereby we introduce the new fermionic operators cl^^/cx^/j, obeying the 

Cu ' anticommutators ({a, b} = ab + ba) 



o 



Then, we write the usual c\./cx as the superposition 



A'o 1 Na 



1^- \ The usual anticommutation relations of c\./ Cx are preserved, while 

l>, ^ 

f— ^ ' {Ca;,Cj/,i/} = , {Cx,Cy^,} — —r= . (3) 

We also introduce 

KN ' ^i,k,i7.S — '^i.k.ij + Wi,k,a C^^k,cr,S '^i.-k.a,S ~^'^+.i,k,a C^^k,a,S '^j.-k,a,S ~^'^-,i,k,<y C^_i.^^_g Cj,.^^g . (4) 

■ - - BJ 1^ g is a, bosonic operator, creating spin triplet pairs of fermions (for singlet pairs c.f. (I)), and {i,j) = {(1, 2), (2, 1)}. 

Henceforth we divide the momentum space into two parts, say k > (sgn(fc) — +) and fc < (sgn(fc) = — ). For 
fc > we form the following multiplet of Bj ^ ^ 's 

^^k = ^l,fc,t,'5=l ^l,fc,4.,5=l ^2,k,t,S=2 ^2,k,S=2 ■ (^) 

This multiplet creates all states with momenta ±fc, and we take No — 2 (c.f. (I) also). 

The new index allows for the bookkeeping of a superposition of states of a given particle, i.e. same x = {i, k,a}, 
without the difficulties due to entanglement within the multiplet Mj., if the index were removed. In that case, the 
treatment of the coherence factors u^w,! is prohibitively complicated. 

This index is "ghostly", i.e. it appears only within the state \'i') belov^r. It does not appear in any 
observable quantity! (C.f. the Faddeev-Popov ghosts in gauge field theories, where the ghosts disappear from 



the final physical outcome.) We note that there is no change whatsoever implied in the Hamiltonian or in the 
representation of any observable. 

Now we introduce the dis-entangled state 



k>Q 

Note that all Blag's in |^) commute with each other. 

|\E') generalizes |\E'bcs) and sustains superfluidity. This wavefunction makes sense for two or more fermion species, 
with an interaction between different species. It can obviously be generalized for three or more fermion species. 
Moreover, a similar wavefunction using the new quantum index can be written in the real space representation 
instead of the momentum space one. At the moment, it is not possible to say how close |^) is to experimental reality. 
It does represent a very promising avenue though, as can be seen from the discussion which follows. 

It allows for inequivalence between spin up and down fermions. Plus, it allows for the "exact" variational treatment 
of a wider class of Hamiltonians than sheer BCS type, e.g. comprising interaction and hybridization between different 
fermion species, in the well known manner of the BCS-Gorkov theory [3|,y|. 

The normalization (^|^) = 1 implies 

\Ui,kA'^ + \Wi^k,a\^ + |t+,i,fe,aP + \t-,^,kA^ = 1 . (7) 

Fermion statistics yields Wi,^k.a = ^Wi^k,a- 

Two fermion species at zero temperature. Calculations are straightforward for the matrix elements derived 
from l^t). E.g. for 2 fermion species with dispersions ei^k,cr = ^i.-k.a and for fc > we have 

\/^ci,fc,t Ml |0) = (wi,fc,o- 4,~k,i,s=i+'t+s.k.t 4,-k,i,5=i) ^l,-fc,fci bIm.s=2 Blk.us=2 |0) 

Then 

(0| Mk cl,fe,^ci,fe,t Ml |0) = ^(ki.fe,tl' + \t+,iM? + \t-,2.kA?) ■ (9) 

(0| Mk cluA^i,kA Ml |0) = -i-(r , ,^^ w,,kA + t^a.k^^ wl^^{) . (10) 

(0| Mfc C2,-fc,iCi,fc,t Ml |0) = ^K.fe.. ^-,i.fc.- - ulk,a t+,2.k,.) ■ (11) 

Using the commutativity of Bl^. ^ ^'s and generalizing the previous equations, we obtain ((C) = (^|C|^)) 

,a^{cU,C,^k,a)^^{\w^,k,a? + \t~l.,^.kA^ + \tl.,3.kA^)) > (^, j) = (1, 2), (2, 1) , (12) 



Likewise, 



and 



TT'i^k^u \~r_k..rr~i-,'i,'J I j\j 



CkM = {4,k,a^3,k,a) = -7r(Kfc,a th,^,k,a + Wj^k,a i*l^,j,kA ' (1^) 



1 



lk,a = {Cj-k,aCt,k,a) = -rp«,fe,^ *ifc,i,/c,<T " Uj,fe,^ t-i^ j.fc.a) , (14) 



with Ik — — sgn(fc). 

A general Hamiltonian for two fermion species interacting via intra-species potentials Vi^2 and via an inter-species 
potential Fq^ and hybridizing via /ifc, is 

H =Y1 ^»^fc.'^ 4,k,a(^^,k,o +^hk (cI,fc,^C2,fc,a + 4,fc,^Ci,fe,^ j (15) 

i.k.a k,{T 

i,k,p,q,(T,(T' k,p,q,a.a' 



with z = 1, 2, £,i^k,cT = U.k.a — fJ-i.tr and fii,a the chemical potentiaL 
Considering * above, we have for (H) = {^\H\^>), ao = I /No 



w = E ^^ 



1 






i,k,a ni^k,a + 2_^ ^k (Cfe.CT + Ck.a) ^ 9 
i,k,a k^a i,k,p,(T 



r, /^ Vi^k-p ni^k,cr ni^p,„ + ^ E ^*:9=o '^i,fc,cr n-i,p,CT (16) 



2,fc,p,cr 









i,k,p.(7 



/ ^ ^ k—p \k,a Sp^a ~^ -^g— / ^ ^l,k,a ^2,p,(T ~r ^ ^ - «:— p /«;,(t /p^, 
k,p,(T k^p.{7 k,p,cr 



fk-p 7fc,CT 7p,cr ) 



with {i,j) = {(1, 2), (2, 1)}. The first term in the second line is exactly the usual BCS-like pairing term, and the last 
term is the equivalent inter-species pairing term due to Fq. 

The minimization procedure for (iJ) and the finite temperature extension proceed as shown in (I). 

Below we present numerical solutions of eqs. (19)-(21) of (I), i.e. the ground state found assuming a real |$) 
and, of course, real superconducting gaps A^ j,. We considered 2 species of electrons in a 2-d lattice with dispersions 
ti^k = —2ti{coskx + cosky) — 4:t^ cos kx cos ky, k^^ky = [— 7r,7r]. We took i,; ^l,t[ = —0.35 and hybridization hk = 0. 
rii is the filling factor of the respective electrons. For the intra-species potential we consider the realistic non-separable 
form (and thus harder computationally) 



V,{q) = V, sin2(q,/2) siv?{qy/2) , 



(17) 



which is peaked at Q = (±7r,±7r). For the inter-species potential we consider two different forms, namely (A) 
Fq = Vb[cos(q:r/2) + cos(gy/2)] and (B) Fq = Vb sin^(ga;/2) sin^((7y/2), i.e. of the same type as the intra-species 
potential. E = (H). An 80x80 discretization of the Brillouin zone was used. 

Vi = l/2=5, ^0 = 1 



m =0.91,n2 =0.809 
Hi = n2 = 0.91 



E = -0.35183 , A?,;, : ^(Ca), A^,^ : ^(Ca) 
E = -0.40270 , A;_fc : D 



E = -0.35660 , A'lk ■ D, A^^^ ; D 
E = -0.40801 , NS 



Vi = l/2 = 5, T/o = 5 



m =0.91,n2 =0.809 
ni = n2 = 0.91 



E = -0.19182 , A?,, : ^(Ca), A°,,. : S{C2) 
E = -0.23120 , A°, : S{±) 



E = -0.35656 , A?,,. : S, A^,, : S 
E = -0.40801 , NS 







Vx = V2 = 10, Fo = 


-- 1 




ni =0.91,n2 =0.809 


E = 


-0.37624 , A;,fc : S, A^,^ : ^'(Ca) 


E = 


-0.38283 , A?,fe : S, A^^ : D + 5(6*2) 


ni = n2 = 0.91 




E = -0.43078 , NS 




E = -0.43726 , NS 


ni = nz = 0.809 




E = -0.32175 , A^,fc : D 




E = -0.32663 , A°fc : S 



Vi = V2 = 10, Fo = 5 



ni =0.91,n2 = 0.809 


E-- 


= -0.17855 , A?,, : S(±), A«,, : S 


E = 


-0.38226 , A?,fc : S, A",^ 


D{C2) 


ni — n2 = 0.91 




E = -0.29048 , A°fe : S(±) 




E = -0.43726 , NS 




ni = n2 = 0.809 




£=-0.13141 , A°fc :S(±) 




E = -0.32718 , NS 





TABLE I: Ground state as a function of the parameters shown. The 2nd column corresponds to Fq = 'K)[cos(gi/2) -|-cos(g;,/2)] 
(A) and the 3rd column to Fq = Vo sin^(ga;/2) sin^(gy/2) (B). S : extended s-wave gap, D : d^2_„2-wave gap, S + D : 
combination of S & D, NS : normal state, A^j. = 0, C2 : symmetry of the gap (C4 otherwise), ± : both signs of the gap are 
present. 
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